Crucial Role of Quantum Entanglement in Bulk Properties of Solids 
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We demonstrate that the magnetic susceptibility of strongly alternating antiferromagnetic spin- 
1/2 chains is an entanglement witness. Specifically, magnetic susceptibility of copper nitrate (CN) 
measured in 1963 Q cannot be described without presence of entanglement. A detailed analysis 
of the spin correlations in CN as obtained from neutron scattering experiments provides micro- 
scopic support for this interpretation. We present a quantitative analysis resulting in the critical 
temperature of 5K in both, completely independent, experiments below which entanglement exists. 



Entangled quantum systems can exhibit correlations 
that cannot be explained on the basis of classical laws. 
Since the birth of quantum theory, such correlations 
have been used to highlight a number of counter-intuitive 
phenomena, such as the Einstcin-Podolsky-Rosen para- 
dox |3j or quantum non-locality |4j - the conflict between 
quantum mechanics and local realism as quantified by 
violation of Bell's inequalities. In recent years entangle- 
ment was realized to be a crucial resource that allows for 
powerful new communication and computational tasks 
that are not possible classically [5j. 

The existence of quantum entanglement is generally 
not seen beyond the atomic scale. Only very recently en- 
tanglement experiments were realized with increasingly 
complex objects, either by entangling more and more sys- 
tems with each other 0, 0, 0, , or by entan glin g systems 
with a larger number of degrees of freedom [lfj . Moving 
towards demonstration of entanglement at even larger 
scales will tackle the question on limits on mass, size 
and complexity of systems that still can contain entan- 
glement. The usual arguments against seeing entangle- 
ment on macroscopic scales is that large objects contain 
large number of degrees of freedom that can interact with 
environment thus inducing decoherence that ultimately 
lead to a quantum-to-classical transition. 

Remarkably, macroscopic entanglement in solids, that 
is, entanglement in the thermodynamical limit of infi- 
nite large number of constituents of solids, was theoreti- 
cally predicted to exist even at moderately high temper- 
atures 0, 0, 0, Q, 0] . Recently, it was demonstrated 
that entanglement can even affect macroscopic thermo- 
dynamical properties of solids 0, 0, 0, 0, 0, Eof . 
such as its magnetic susceptibility or heat capacity, al- 
beit at very low temperature (few milikelvin) and only for 
a special material system - the insulating magnetic salt 
LiHo x Yi_ x F4 [2jj. This extraordinary result shows that 
entanglement can have significant macroscopic effects. 

Nevertheless, it is an open question whether or not 
there are macroscopic thermodynamical quantities that 



are entanglement witnesses for broader classes of solid 
state systems. Entanglement witnesses are observables 
which have positive expectation values for separable 
states and negative ones for some entangled states [22|. 
Finally, a direct experimental demonstration of entan- 
glement, in the form of a determination of correlations 
between microscopic constituents of solids, remains an 
experimental challenge. 

Here we show that bulk properties of solids can be en- 
tanglement witnesses. We use already published experi- 
mental results of both microscopic structure and macro- 
scopic properties of the spin-1/2 alternating bond anti- 
ferromagnet CN (Cu(N03)22.5D20), to demonstrate the 
existence of entanglement in this type of solid. In the 
first approach we analyse experimental results of neutron 
scattering measurement of CN obtained in 2000 and 
show that they provide, for the first time, a direct ex- 
perimental demonstration of macroscopic entanglement 
in solids. The experimental characterization of the dy- 
namic spin correlations for next neighbouring sites en- 
ables us to determine concurrence |2ll - a measure of 
bipartite entanglement - and show the existence of en- 
tanglement at moderately high temperatures (as high as 
5 Kelvin). In the second, parallel, approach we show 
that magnetic susceptibility at zero magnetic field is a 
macroscopic thermodynamical entanglement witness for 
the class of solid states systems that can be modeled by 
strongly alternating spin-1/2 antiferromagnet chain. We 
then show that the measured values for magnetic suscep- 
tibility of CN in 1963 [lj imply presence of entanglement 
in the same temperature range (below 5 Kelvin). 

CN is an accurate realization of strongly alternat- 
ing one-dimensional antiferromagnetic Heisenberg spin 
chain. The corresponding spin Hamiltonian is given by 

&2j ' &2j+l + J2&2j+1 ■ S2J+2), (1) 

3 

representing pairs of spins that are alternately coupled 
by strong intradimer 3\ « 0.44 mcV and weak inter- 



dimer J 2 « 0.11 meV coupling 23]. This can be de- 
scribed by a model of antiferromagnetically coupled spin 
pairs - dimers - which are themselves coupled by weaker 
antifcrromagnetic interaction. The model has a highly 
entangled nontrivial spin-0 ground state |24l Ejf and for 
all < J2 1 J\ < 1 has a gap of the order of A » J\ to 
the first excitation, which is a band of spin-1 excitations 
(magnons). However, because here J2I 'Ji~0.24 @,H3 , 
it is useful to think of CN as a chain of uncoupled spin 
dimers. Each dimer then has a singlet ground state and 
the triplets are the degenerate excited states. The exis- 
tence of the energy gap gives the first estimate for per- 
sistence of entanglement for temperature range below 
T« Ji/k~5K, where fe is the Boltzmann constant. 

Next we describe the main experimental results of 
Ref. 0- We will follow the discussion given there. 
CN has a monoclinic crystal structure with space group 
/12/cl and with low-temperature parameters a = 16.1 A, 
b = 4.9l, c = 15.81, and j3 = 92.9°. The vec- 
tor connecting dimers center to center is Uo = [lll]/2 
for half the chain, and u' = [III]/ 2 for the other 
half. The corresponding intradimer vectors are di = 
[0.252, ±0.027, 0.228], respectively. In the experiment the 
neutron scattering intensity was measured in the temper- 
ature range 0.3LK" < T < 7.66K (i.e. O.O6J1 < kT < 
1.5Ji) as a function of energy transfer huj and wave 
vector transfer Q. Count rates were normalized to in- 
coherent elastic scattering from the sample to provide 
absolute intensity I(Q,u) = || F(Q)| 2 2S(Q, uj). Here 

.9 = V(9b + 9±)/ 2 = 2 - 22 with 9b = 2 - 31 and 9± = 2.11 
that show small anisotropy of <?-factor along and perpen- 
dicular to the crystallographic direction b [2^], F(Q) is 
the magnetic form factor for Cu 2+ j^], and S(Q,us) is 
the scattering function |27|. 

The direct link between the microscopic structure as 
given by the correlation function between spins and the 
intensity of inelastic neutron scattering is given by an 
exact sum rule (the first cj-moment of scattering cross 
section) |28| 

/+00 
i5(Q, u>)du) 
-OO 

= -i^Jd(S -S d )(l-cosQ.d), (2) 
6 d 

where {d} is the set of all bond vectors connecting a spin 
to its neighbours, S(Q,uj) is single site normalized and 
(S • S d ) = (S$S%) + (S y Q S y d ) + (S^S z d ) is the sum over 
correlations for three orthogonal directions x, y and z. 

The intradimer correlation (So • S dl ) between next 
neighbouring spins was extracted from the global fits of 
the following phenomenological form for S(Q,uj) to the 
complete data set at each of temperatures (it was used 
more than 1000 data points per parameter): 
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1: The temperature dependence of spin correlation func- 
tion (So - Sd, } for neighbouring sites in CN. The figure is taken 
from Ref. J. The open circles correspond to the tempera- 
tures at which measurements were performed, the solid black 
line were obtained from fits in Ref. 0- The solid red line 
is from this work. It distinguishes temperature ranges with 
and without entanglement in CN. The critical temperature is 
around TJf xp ~ 5K and 0.25 on the y-axis is the maximal value 
for (So • Sdj) achievable with separable states. For T <T^ xp 
concurrence is given by C = — 2(So ■ Sdx) ~ (1/2) and has the 
same functional dependence on temperature as given in the 
figure. For T>T* xp , C vanishes. 



Here /3 = l/(fcT), H{u)q stands for expression J2J, f(E) 
is a normalized spectral function and e(Q) is the disper- 
sion relation. Equation @ represents the "single mode 
approximation" which is valid for sufficiently low tem- 
peratures |2j| . The dispersion relation is used in the 
variational form based on the first order perturbation 
theory |H 



e(Q) = Ji - ~ ^2 J u cos Q • u, 



(4) 



S(Q,w) 



1 



e(Q) l-exp[-j9e(Q)] 



f[Hw-e(Q)]. (3) 



where the vectors {u} connect neighbouring dimers cen- 
ter to center both within and between the chains. 

At T = 0.3K they used f(E) = 5(E) and obtained 
the global fit with an overall factor and four exchange 
coupling constants (J\, J2 and the constants Jl and Ju 
for coupling between the chains 0) in Eq. as the 
only fit parameters. The fit gives (So • S dl ) = —0.9(2) 
in agreement with the value of -3/4 for dimers in the 
singlet state. To perform the fit at higher temperatures 
the spectral function was replaced by a normalized Gaus- 
sian with half-width-half-maximum T(q) = Tq + 5^- cos q, 
where q = Q • u is wave vector transfer along the chain. 
Also, the dispersion relation 10} was replaced by the form 
e (Q) = Ji— n {T)\ J2u "^u cos Q' u j where renormalization 
factor n(T) was introduced to account for finite temper- 
atures. The prefactor for global fits at each temperature 
yields the intradimer spin correlations (So • S dl ) as given 
in Fig. ^ The correlations decreases with temperature 
increase due to mixing of the triplets with the singlet 
state of each dimer. 

We will now show that the values estimated for the cor- 
relation function can only be explained if entanglement 
is present in the solid. We first show that | (So • S dl )| is 
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an entanglement witness. The proof is based on the fact 
that for any product state of a pair of the spins one has 

|(s -s dl )| = |<^)(SS 1 > + (S2[)(S3 1 ) + <S5)<55 1 )I 

< |S ||S dl |<l/4. (5) 

Here and throughout the paper we choose units to be 
consistent with those of Ref . |2| ; spin is expressed in units 
of 1/2; H = 1. The upper bound was found by using 
the Cauchy-Schwarz inequality and knowing that for any 
state |S| = ^(S x ) 2 + (Sv) 2 + (S z ) 2 < 1/2. The proof 
is also valid for any convex sum of product states of two 
spins (separable states): p = £ fc w kP\®Pk with £ fc Wk = 
1. In Fig. the value of 1/4 for — (So • S dl ) distinguishes 
temperatures ranges with and without entanglement in 
CN. The critical temperature is found to be T^ xp « 5K. 
The reported error in the correlation function (A = 0.2 
at T = 0.3K 0) implies an error of around AT » IK in 
the critical temperature. 

Within the model of uncoupled dimers the isotropy 
of Heisenberg interaction in spin space ensure that 
(SqS^) = (SqS^ ) — (SqS^ ) and concurrence is given 
by C = 2max[0, -(S • S dl > - (1/4)]. Similarly Bell's pa- 
rameter - the quantum value of Bell's expression in the 
Clauser-Horne-Shymony-Holt inequality |3lj - is given by 
(8\/2/3)|(So • S dl )|. It is higher than the local realistic 
limit 2 at temperatures below T^ xp . Apart from rescaling 
both the temperature dependency of concurrence and of 
Bell's parameter have the same functional form as the 
correlation function in Fig. ^ 

The temperature dependence of (So ■ S dl ) in Fig. 
is in a good agreement with a model of uncoupled 
dimers for which (S • S dl ) = — (3/4)An.(/3Ji), where 
An(J3Ji) = (1 - e- /3Jl )/(l + Se-f 3 -' 1 ) is the singlet-triplet 
population difference Q- Within the model the theoret- 
ical temperature dependence of concurrence is given by 
C = max[0, (1 - 3e _/3 Jl )/ ( 1 + ie' 13 Jl )] , as first obtained in 
Ref. [Ill \l2\. The theoretical value for the critical tem- 
perature Tf = Ji/(fcln3) = A.6K is in excellent agree- 
ment with the value estimated from the experiment. 

We now proceed with our second approach. We will 
analyse experimental results of a magnetic susceptibility 
measurement of CN to show that the values measured 
at low temperatures cannot be explained without entan- 
glement being present in the system. This will be based 
on a general proof that magnetic susceptibility of any 
strongly alternating antiferromagnetic spin-1/2 chain is 
a macroscopic thermodynamical entanglement witness. 

When the system is in thermal equilibrium under a 
certain temperature T, its state is p = e~ H ' kT /Z, where 
Z = Tr(e~ H ' kT ) is the partition function and H is the 
Hamiltonian. The magnetic susceptibility along direction 
a is defined as Xa = (d(M a )/dB) = (g 2 p 2 B /kT)(((M a ) 2 ) - 
(M a ) 2 ), where (M a ) is magnetization along a, M a = 
Y^jSj 1 , B is external magnetic field, g is g- factor and 
Pb is the Bohr magneton. Because of the isotropy of 



the Hamiltonian in spin space magnetization at zero- 
field vanishes for any temperature. This implies the 
following form for magnetic susceptibility at zero-field: 
X a = {g 2 p 2 B /kT){{M a ) 2 ) B=0 = (g 2 p 2 B /kT) £. . {S?S?) « 
(g 2 p 2 B N/kT) [(1/4) + {SgSgJ]. Here we assume that cor- 
relations between all spins that are not nearest neighbors 
are negligible compared with (SqS^ )■ This approxi- 
mation is valid for our case of strongly alternating spin 
chains at low temperatures. It is important to note that 
apart from the weak anisotropy in the g factor, the mag- 
netic susceptibility at zero-field is isotropic, i.e. \x — Xy — 
Xz = X- Thus, if we sum the values of magnetic suscepti- 
bilities over the three orthogonal directions x, y and z in 
space we obtain x ={g 2 P 2 B N/kT)[{l/A) + ((S • S dl )/3)], 
where the mean value (...) is taken over the thermal state 
at B = 0. However, because one has | (So ■ S dl )| < 1/4 for 
any separable state (Eq. J5J), the magnetic susceptibility 
for such states is limited as given by 



X> 



kT 6' 



(6) 



Thus, a violation of this inequality necessarily detects 
entanglement in the system. 

In Ref. 0] magnetic susceptibility of CN was measured 
on single crystal in 0.4 — 4.2K and 14 — 20K ranges of 
temperature. The measurement method was based on 
a mutual inductance bridge working at 275 Hz. The 
susceptibility was found to have a rounded maximum 
at 3.2K dropping very rapidly below this temperature 
approaching zero at vanishing temperatures, as given in 
Fig. |21 Such behaviour is typical for alternating spin 
chains. Thermodynamical entanglement witness is 
represented by the red solid line in Fig. [3 The measured 
values of magnetic susceptibility below the intersection 
point of the curve representing the witness and the ex- 
perimental curve cannot be described without entangle- 
ment. The critical temperature is around T^ xp w 5K. 
This is in excellent agreement with the value estimated 
from the neutron scattering experiment in spite the fact 
that different samples were used (authors of Ref. 1] noted 
possible variations of the physical state of the sample due 
to high hygroscopy of CN), the two experimental meth- 
ods test entirely different physical quantities and are al- 
most 40 years apart (we note that no experimental error 
of magnetic susceptibility measurement was reported in 
Ref. U). 

In conclusion, we show that neutron scattering experi- 
ments YM directly and measurement of magnetic suscepti- 
bility [l| indirectly demonstrate presence of macroscopic 
quantum entanglement in a solid (cooper nitrate). We 
believe our results indicate that entanglement may play 
a broad generic role in macroscopic phenomena. 

The question of having macroscopic entanglement is 
not only fascinating in its own right but it also has a fun- 
damental significance as it pushes the realm of quantum 
physics well into the macroscopic world, thus opening the 
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FIG. 2: The temperature dependence of magnetic suscepti- 
bility of powder CN (triangles) and a single-crystal CN mea- 
sured at low field parallel (open squares) and perpendicular 
(open circles, crosses, filled circles) to the monoclinic 6 axis. 
The data and the figure are from Ref. Q. The solid black 
curve is the theoretical curve for a dimer re-scaled for the 
amount of noise estimated from the experiment. This noise is 
computed as the ratio of the maximal experimental value (av- 
eraged over crystal data) and the maximal theoretical value. 
The solid red curve represents the macroscopic entanglement 
witness © and is re-scaled in exactly the same way. The 
intersection point of this curve and the experimental one de- 
fines the temperature range (left from the intersection point) 
with entanglement in CN. The critical temperature is around 
rpemp ^ Note that the entanglement witness will cut the 
experimental curve (and hence yield a critical temperature) 
independently of a particular re-scaling procedure. 



possibility to test quantum theory versus alternative the- 
ories well beyond the scales on which theirs predictions 
coincide. It also has important practical implications for 
implementation of quantum information processing. If 
the future quantum computer is supposed to reach the 
stage of wide commercial application, it is likely that it 
should share the same feature as the current (classical) 
information technology and be based on solid state sys- 
tems. It will thus be important to derive the critical 
values of physical parameters (e.g. the high-temperature 
limit) above which one cannot harness quantum entan- 
glement in solids as a resource for quantum information 
processing. 

We consider our work to imply that many experiments 
performed in the past may still hide new and interesting 
physics. Experiments of Berger et al. 1] from 1963 and 
Xu et al. |2| from 2000 used here are exemplary. Inter- 
estingly, the 1963 experiment was performed long before 
any serious attempts to measure entanglement begun in 
the seventies. It is not even impossible that we are able 
to find a result from which we can infer the existence of 
entanglement and which appeared well before this con- 
cept was conceived by Schrodinger in 1935. 
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